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Traditional approaches to designing quantum memories presume a complex classical 
apparatus operating in parallel with the memory to detect and correct errors as they 
arise [U [2]. Typically, the mechanics of this essential half of the memory system are 
hardly considered at all, taking a back seat in the analysis to the qubits in which the 
information is stored. This assumption of a perfect, classical overseer glosses over the 
inherent technological mismatch of typically fast, nanoscale and cold quantum systems 
with slow, mesoscopic, and/or hot classical ones. The extreme performance demands 
for quantum information processing should motivate solutions that weigh classical and 
quantum resources wholistically and utilize system models as close to physical mechanics 
as possible. 

In [3] we described an approach to designing protected quantum memory devices 
naturally suited to nanophotonic implementations that is technologically homogeneous 
and requires no "off-chip" oversight, simply "power" in the form of cw laser inputs. The 
particular design example in [3] emulates the bit- or phase-flip quantum error correcting 
(QEC) code [U |2] with a single logical qubit stored in the collective state of three multi- 
level "atoms," each strongly coupled to three different optical resonators. Two more 
cavity quantum electrodynamical (cQED) systems serve as the binary "controllers" 
[1] for these qubits. These five cQED devices are connected by a network of single- 
mode waveguides and beamsplitters that when appropriately powered by cw laser inputs 
continually and simultaneously cause the internal states of the controllers to reflect the 
joint parities of the qubits and corrective feedback on the qubits to be implemented 
according to the controller states. Essentially a bit-/phase-flip QEC code, the design in 
[3] can only protect the logical qubit against either single-qubit bit-flip or single-qubit 
phase-flip errors, depending on the configuration. To emphasize how straightforwardly 
the approach scales to more powerful QEC codes, in this article we describe, model 
and simulate an autonomous nanophotonic network that emulates the 9 qubit Bacon- 
Shor subsystem code [5j E], the smallest of a class of naturally fault-tolerant QEC 
codes capable of protecting a single logical qubit from arbitrary single-qubit errors (see 



Appendix A for a description of this code). 



While these networks are conceived of intuitively, our analysis arises from a 
mathematical framework of open quantum optical systems that takes quantum field 
theory as its physical basis [7J |8] but most closely resembles a quantum generalization 
of electrical circuit theory [9J UP] . Although rigorous stochastic dynamical modeling 
is often unfamiliar to physicists, this formalism [TT] is very physically intuitive once 
internalized and sufficiently flexible to describe the continuous-time dynamics of most 
systems foreseeable in quantum optical networks. In our case, cQED devices are modeled 
as distinct Hamiltonian systems that couple weakly to free, bosonic fields that scatter 
off of each device in series and in parallel along paths set by the single-mode waveguides 
linking the network devices [9]. As our network operates "autonomously," without 
any user monitoring necessary, the equation of motion that describes the closed loop 
time-evolution of the devices (i.e. the dynamics after tracing over the field degrees of 
freedom) is a deterministic master equation [8j |3]. Moreover, as all the Hamiltonian 
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Figure 1. a) Schematic of nanophotonic network capable of implementing the 9 
qubit Bacon-Shor QEC code. CW coherent field inputs that probe the "Z" and "X" 
syndromes of the memory qubits, Qij, enter from the middle of the bottom and 
left-hand side, in blue and green, respectively. After traversing the memory qubits, 
the phases of these fields represent measurements of the four syndrome generators. 
Through interference with four more cw "local oscillator" laser inputs on beamsplitters 
and interaction with four "relay controller" qubits, Ri, these phases effectively control 
the relays' internal states. The relay internal states then direct four "feedback" 
cw inputs towards the memory qubits. When two red (orange) feedback beams 
simultaneously illuminate a memory qubit, coherent Pauli-X (-Z) rotations occur 
until a "no-error" syndrome state is recovered, at which point the corrective feedback 
dynamics automatically shut off. b) & c) Example memory and relay cQED input- 
output, internal level structure, and coupled atomic transition schematics, adapted 
from [3]. 

and device-waveguide couplings are constant in time, the network is stationary. Much 
like an electronic operational amplifier with a feedback impedance network, together 
the cQED memory and controllers represent an integrated, self-stabilizing system that 
simply requires DC "power" to function. 

As in [3] , the memory storage qubits are physically realized by multi- level "atoms" 
with two ground states that represent the spin-up and -down states of an ideal qubit. 
When an excited state couples to only one ground state (in some basis) via an electric 
dipole transition that is degenerate with and strongly coupled to a mode of a single- 
sided optical resonator, then, in appropriate limits, an on-resonance cw laser beam may 
scatter off the resonator without dissipation or perturbing the qubit state, but will 
acquire a 7r phase shift upon reflection if the atom is in its coupled ground state, or 
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no phase shift if in its uncoupled state. This probe interaction approaches this ideal as 
the atom-mode and waveguide-mode coupling rates dominate all other dynamical rates 
(dissipative process rates, in particular) [12], which is the natural limit to consider in 
a "small-volume," quantum nanophotonic context. The identification of the physical, 
multi-level system with the ideal, dissipation-less system is formalized by applying these 
limits and an adiabatic elimination theorem built to work in our quantum optical 
framework [13] . Whereas in the bit-/phase-flip network [3], only one type of probe 
interaction was required per qubit, in this more powerful network, each memory qubit 
is equipped with the necessary atomic transitions and single-sided resonator modes 
required to probe the atomic ground state in a Pauli-Z basis with one external probe 
field and in a Pauli-X basis in another. Similarly, the controller cQED devices approach 
their ideal functionality in "small-volume" limits equivalent to those of the memory 
devices In these limits, each operates as an optical set-reset relay that directs the 
optical power input of one channel out of one of two output channels depending on the 
relay's internal state. This binary internal state is in turn controlled by the presence or 
absence of optical power entering two additional optical inputs. 

Using these components and the intuition gained in [3j, it's in fact immediately 
clear how one could design a self-correcting nanophotonic network that emulates the 
9 qubit Bacon-Shor QEC code [5]. In addition, one may also immediately write 
down the full master equation that governs the internal dynamics of all memory and 
controller devices of this network without having to resort to a quantum circuit model, 
suggesting the general applicability of the approach (a laborious, but straightforward 
network calculation confirms the master equation; see below). In figure [1] we sketch 
the nanophotonic design. Two probe laser inputs first travel down the middle row and 
column of a 3 x 3 grid of memory qubits. At the end of the middle row or column, both 
probes are split on a beamsplitter with each of the now four probe beams traveling back 
along the remaining rows and columns, respectively. The frequency and polarization 
of the laser probe traversing the rows (columns) is such that upon reflection from each 
memory cQED device, the probe picks up either a tt or phase shift according to the 
atomic ground state in a Pauli-X (Pauli-Z) basis. Thus after traversing the memory 
qubits, the phases of the four probe lasers become entangled with the memory qubits 
via controlled-NOT interactions (in a qubit Z- or X-basis) - essentially as four ancilla 
qubits would be to extract the error syndrome in a measurement-based QEC scheme 
[U El |6] . By interfering with four local oscillator lasers, the phases of the four probes 
set the internal state of four relays, which thus play the role of the classical register 
in a measurement-based scheme; the four relay states together represent the quantum 
memory's error syndrome PQ [2]. As in [3], these relays control the routing of four more 
input "feedback" lasers that perform corrective, unitary rotations on the ground state 
of the memory qubits via Raman interactions: when a single qubit is simultaneously 
illuminated by feedback beams emitted by both "Z-syndrome" ( "X-syndrome" ) relays, it 
undergoes coherent Pauli-X (Pauli-Z) rotations between its ground states (for simplicity 
in this much larger network, we assume Stark-shift compensation mechanisms are in 
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Figure 2. Schematic of the "zig-zag" configuration of the Z-probe network only. Due 
to the subsystem structure of the code, dotted connections are capable of suffering 
calibrated waveguide loss without disruption to the logical information stored in the 
memory qubits (see | Appendix A[ ). The "zig-zag" configuration for the X-probe network 
follows analogously, essentially rotated by 90°. 

place [3]) until the memory qubits - followed by the controllers - recover their "no- 
error" state, and automatically shut off the corrective feedback. Although errors may 
occur to any of memory qubits, corrective feedback is applied to only a small subset of 
the memory qubits, as in figure [U due to the subsystem structure of the code fJU [5] (see 



Appendix A). Moreover, as the internal states of the relays are continually reenforced 
by the cw probes, errors in the relays' atomic states should be self-corrected by the 
network. 

The subsystem structure of the Bacon-Shor code [HI E] may be further leveraged 
in a slightly modified network configuration. Although doing so would require 
a significantly more intricate waveguide network and/or cQED devices, instead of 
threading the probe fields directly along rows and columns, a zig-zag configuration, 
as in figure [2J should aid network robustness against optical waveguide loss. Whereas 
photon loss almost anywhere in error-probe waveguide depicted in figure Q] effectively 
causes network- induced "errors" (which sometimes can be ameliorated by the correction 
network), logical fidelity in the "zig-zag" network should be immune to any waveguide 



loss in half of the probe network connections [6] (see Appendix A ). We will consider 
these losses and other critical robustness concerns in depth a later publication, however, 
and so for our purposes here, both networks represented in figures[T]and[2]have equivalent 
dynamics. 
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The quantum circuit model approach to deriving the master equation dynamics 
of the network depicted in figure [I] proceeds by appropriately "connecting" individual 
quantum optical device models in series and parallel along freely-propagating optical 
field modes. In our quantum stochastic differential equation (QSDE) formalism [11], 
each device is characterized by a triple (S,L,H), where H is the Hamiltonian of the 
device's internal degrees of freedom; L is a vector of coupling operators, each element 
of which causes these internal variables interact with a particular free field mode; and 
S is a scattering matrix for the input and output field modes, whose coefficients can in 
general be operators on the internal state space of the device [9]. Composite devices 
may be constructed by feeding the field output (s) of one device directly into another 
using the "series product" operation, <: (S 2 i, L 2 i, H 2 i) = (S 2 , L 2 , H 2 ) < (Si, L±, Hi) 
where S 2 i = S 2 Si, L 2 i = L 2 + S 2 Li, and H 2 i = Hi + H 2 +\m[L\S 2 Li\, representing 
a system in which the output (s) of device 1 feed(s) into 2. Devices may also be 
connected in parallel using the "concatenation product," EH, by which, for example, two 
single-input /-output devices may be aggregated into one dual- input /-output device: 
(diag{S 2 ,Si},{L 2 , L x }, H 2 + H x ) = (S 2 ,L 2 ,H 2 ) EH (Si,Li,Hi). The catalog of the 
quantum devices in our network (each with an associated triple) includes [3] the n- 
field identity element, J n ; cw laser inputs with amplitude 7, W 7 ; 50/50 beamsplitters, 
B; the Z-probe (X-probe) aspect of each memory qubit, Qf, (Qfj)] the set-reset control 
of the relays, R\\ the feedback routing control by the relays, R\\ and the individual 
drive-branches of each Raman-X (-Z) interaction in of the memory qubits, Qf^ (Qfj) 
k G TL 2 . With these in hand, the appropriate network model may be essentially read off 
from diagrams like figure [H For example, the entire, closed loop QSDE model for our 9 
qubit Bacon-Shor network may be constructed using these individual elements by 

G z = (R\ EH R c 2 ) <(BSB)< ((Q* < Q z %x < Qfj EH h EH (Q^ 3 < Q% >3 < Qfj) < 

((B EH 2 h) EH h) < (JfiZ 2 < QI2 < QI2 < EH W a EH h EH W°^j 

G n = (Qf t l EH Qfj EH Qfj < (h EH B) < ((i?f < (W? EH /j) EH I x ) 

G f2 = (Q*l EH Qg EH Qfj < (h EH B) < < (W? EH jj) EH I^j 
Gp = (R c 3 EH Rl) < (B EH B) < ((Q* < Q* 2 < Qfj EH h EH (Q* < < Qfj) < 
((B EH 2 h) EH h) < ((q* < Qf 2 < Q* < EH W a ffl h ffl W°^j 

Gf3= (Qll EH i EH Qfj < (h EH B) < ((i^ < (W* EH EH 

G /4 = (Qg EH gg EH Qf*) < (7 X EH £) < (Jr{ < (W p ffl /j) EH 

G network = G p EH G71 ffl G72 EH Cr^ ffl G73 EH G/4 (1) 

where G z (Gp) describes the "Z" ("X") syndrome extraction network, Gfi describes 
the feedback network controlled by relay i, and ffl 2 is a "padding operator" related 
to the concatenation product [3] (necessary for proper field indexing), and G networ k 
represents the (S, L, H) triple of the entire system, describing all but the error dynamics. 
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While the symbolic calculation of the analogous G netW0 rk was done by hand in [3], 
this tedious calculation may be automated using simple symbolic matrix manipulation 
scripts in programs such as Mathematica, freeing a quantum network designer to easily 
edit and adjust systems essentially at the schematic level. As in |3J, the final step 
in the calculation, which yields a relatively low-dimensional total network G ne t WO rk, is 
the application of an algorithmic adiabatic elimination procedure [13] that restricts the 
Raman interaction dynamics to the ground states of the memory qubits only 

Including arbitrary single-qubit dephasing of each memory qubit to G netW0 rk 
through simple concatenation products, the closed loop dynamics of p t , the density 
matrix for all 9 memory and 4 relays qubits (physically representing the entire system's 
dynamics after tracing over the field degrees of freedom), is read off from the entire 
system's (S,L,H) triple [HIE], giving: 

pt = - i [if, Pt] + [ L iPt L j ~ \m L hPt}) , with 

h = ft(v / 2Avnrn+ + x 3 , 2 nrn 2 - v^x^n+n, + 



l 3 = (n 2 of &2 + o~2 e^ 2 ) , l 4 = — (-ct 2 + of &2 - Ef &2 ) , 
Ls = ^ fa ** + ' Le = ^ (- ff 3 + of &2 - n 3 + £* 2 ) , 

L 9 _ 17 = VrXij, Lis-26 = VTZjj, L 2 e-35 = vf^ij i,jeZ 3 , (2) 

where {Xij,Zij} are the Pauli-{X,Z} operators on the memory qubit in row i and 
column j, 11^ is the projector onto state |±) of relay i, af T = |±)(=f| for relay i, and 
{ — Oik,j-> Ei&j, — Of&p O i&€ p } = 1 + {—Z^^Z^j, Z^^Z^j, —Xi tJf Xj^, Xi^Xjj,} with index 
* signifying the operator product of all operators acting on some row or column. T 
is the mean error rate of each type of single-qubit error, a is the coherent amplitude 
of the (eventual) four probe lasers (where \a\ 2 has units photons time -1 ), and Q is 
proportional to the optical power in each feedback beam. Comparison with the master 
equation derived in [3] , reveals that these dynamics for the 9 qubit QEC subsystem code 
appear as an almost trivial expansion of that simple bit-/phase-flip network. 

Despite our idealized description of all cQED systems as simple qubits, Eq. [2] still 
represents complex dynamics on a (sparse) density matrix with dimensions 2 9+4 x 2 9+4 . 
Due to the size of the problem, multi-threaded numerical integration of Eq. [2]was carried 
out on a multicore computational server using BLAS routines with OpenMP. Figure [3J 
depicts the time-evolution of the fidelity of the logical information stored in the system, 
initialized with the logical qubit in the +1 Yj, eigenstate, the four gauge qubits in +1 
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Codeword fidelity decay 




Time IT t] 



Figure 3. Decay of the logical codeword fidelity (see |Appendix A| in time for several 
feedback parameters. In the simulation, single-memory qubit X-, Y- and Z-errors occur 
each at average rate T — .1 per qubit. While the amplitude of the probes, are fixed at 
a = 0/8, the feedback "strength" il is varied from (no feedback) to 200. For large 
feedback strengths, the logical storage fidelity is eventually superior to that of a single, 
bare qubit suffering arbitrary dephasing. 



Zci eigenstates (see Appendix A ), and the four relays each in their "no error" |+) state. 
The bitwise-error rate was fixed at V = .1, the probe strength kept at a = Q/8, and the 
"feedback strength" Q varied from to 200. Initially, all fidelity curves drop steeply, 
representing a network "latency" before both probe and feedback networks begin to 
correct errors. Thereafter, fidelity still drops due to the finite probability that multiple 
errors irrecoverably accumulate before correction, but at a retarded rate that decreases 
with increasing feedback strength. For Q > 100, despite the initially steep fidelity loss, 
the decay is slow enough to eventually achieve superior storage fidelity than a single, 
uncorrectable qubit suffering X-, Y-, and Z-errors at the same bit-wise rate. 

We may begin to break down the dynamics of this particular network by calculating 
its performance in response to different types of errors. In figure HI we plot the fidelity 
decay in time with {T,a,Q} = {.1, 100/8, 100}, but with the memory qubits suffering 
only single-qubit X-, Z-, or Y-errors. Systems suffering only X- or Z-errors require only 
"half" of the feedback network to function well and achieve the best storage fidelity. 
The system responds slightly worse to Y-errors. This is expected as Y-flip errors (i.e. 
simultaneous bit-flip and phase-flip errors on a single qubit) requires the two halves of 
the correction network to work in concert to correct the error. That such coordination 
occurs without too much additional fidelity loss suggests that the apparatus is fairly 
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QEC performance for different error classes 
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Figure 4. Comparison of network storage fidelity between systems suffering different 
types of errors. All solid curves are simulated with {V, a, 57} = {.f , 100/8, 100}, except 
for the magenta curve, which has no feedback, as indicated. The blue curve represents 
the storage fidelity when only X- or only Z-errors occur to each memory qubit. The 
only slightly lower red curve comes from a system suffering only Y-errors, which require 
coordination between all aspects of the network. The green curve represents a system 
suffering all three types of error, each at the same bit-wise rate, and has roughy the 
same fidelity decay as a sum of the decays from each type of error acting alone. The 
dashed black line again represents the logical fidelity of a single, bare qubit suffering 
all three types of errors at mean rate T = .1 each. 

unbiased towards the types of errors it corrects. Similarly, when all three errors are 
acting simultaneously, each at mean rate T = .1 per qubit, the additional fidelity loss 
scales roughly linearly with the total error rate in this strong feedback regime, as one 
would expect for an efficiently operating system. 

In conclusion, we have demonstrated that the concept of autonomous nanophotonic 
quantum circuits is sufficiently flexible to be straightforwardly applied in designing 
robust and naturally fault-tolerant quantum memories that emulate measurement-based 
approaches like the 9 qubit Bacon-Shor QEC code. Modeled with continuous-time 
dynamics, our models are built upon true quantum circuit diagrams like figure [IJ in 
which components represent physical devices (not discrete operations) and connections 
represent physical optical fields (not just time-ordered process flows). Future directions 
in our work will focus on developing freely-available software packages that will allow 
users to design, analyze, and simulate similar networks with configurable connections 
and fully-specifiable Hamiltonian devices truly at this schematic level, with only minimal 
expertise in the underlying formalisms. On-going work directly related to these memory 
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networks include analyzing robustness to mis-tuned parameters and optical losses, 
formulating largely automated quantum memory registry "read-in/-out" procedures 
compatible with these networks, and considering protection schemes beyond simple 
emulation of well-known quantum stabilizer codes. 
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Appendix A. The 9 qubit Bacon-Shor QEC code 

This appendix summarizes the measurement-based 9 qubit Bacon-Shor QEC (5j E] 
subsystem code our autonomous photonic network emulates. As in the main article, 
the 9 qubits are arranged in a 3 x 3 grid. In this configuration, a non-Abelian operator 
group is the set of Pauli operators generated by pairs of X operators acting on adjacent- 
row qubits and pairs of Z operators acting on adjacent-column qubits 

T = (XijXi+ij, Z jA Z jii+1 \i e Z 2 ,j E Z 3 ), (A.l) 

where Xij and Zij act as Pauli operators on the qubit in row i, column j and as the 
identity on all others. A subgroup S C T serves as the stabilizer [TJ [2] for the code and 
is generated by X operators acting on all qubits in two adjacent rows and Z operators 
acting on all qubits in two adjacent columns 

<S — (Aj^Aj+i^, Z*^Z* t i + i \ieZ 2 ). (A.2) 

As S is an Abelian group of observables, the ±1 eigenvalues of the generators of S may 
be used to label subspaces of the "H = (C 2 )® 9 9-qubit Hilbert space 

n = @n s (A.3) 

a 

where a is string of the four binary eigenvalues of the stabilizer generators and each Tig 
represents a 2 5 -dimensional space. Each Tig may be further decomposed into subsystems 

H S = UX a ® Hi (A.4) 

where elements from T (which commute with S) act non-trivially on some Hj, and as 
the identity on Tig, encoding four logical qubits, modulo S. This action leaves the H§ 
as single-qubit subsystems, in one of which our logical qubit codeword is stored. We 
can choose this logical qubit to be the one encoded by the operators X L = Xi^ and 
Zl = Z* } i in the a = {+1, +1, +1, +1} stabilizer subspace. The other four logical qubits 
that live in this subspace's 1-C[ =im subsystem may be specified by four commuting pairs 
of anticommuting operators in T: {Xci, Zci\i € Z 4 }. The logical qubits that live in 
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%Z =im are deemed "gauge" qubits that have no information encoded in them by the 
user. The QEC procedure that protects the single logical qubit codeword in %^ =im 
may perturb the state of these gauge qubits without consequence. 

In the measurement-based QEC code, the encoded codeword may be rendered 
impervious to arbitrary single-qubit rotations of any of the nine qubits, although the 
state of gauge qubits may be depolarized. Measurement of the four stabilizer generators 
"discretizes" any partial qubit rotations, projecting the system state into some error 
space Hg, with a = {+1, +1, +1, +1} being the "no error" space. The results of 
these syndrome measurements localize Pauli-X errors down to the column in which 
they occurred, Pauli-Z down to the row, and Pauli-Y down to the qubit. Although 
the precise location of the error is not known generally, codeword recovery is achieved 
by applying a Pauli-X (-Z) to any qubit in the same column (row) as a detected bit- 
flip (phase-flip) error. For example, if a partial Pauli-Y rotation occurs to qubit Q 2 ,2 
(the middle qubit in the 3x3 grid), syndrome measurement may project the system 
into the a = {— 1, — 1, — 1, — 1} subspace, digitizing and indicating the error. The 
encoded codeword may be recovered by then applying the operator X 3;2 Z 2i i, as the 
net, recovery-error operator, X^ 2 Z 2i \Y 2 ^ commutes with S, Xl and Zl- However, 
-^3,2-^2,1^2,2 £ T and thus the recovery operation disturbs the state of the logical gauge 
qubits. Moreover, the error syndrome need not be acquired by directly measuring the 
6-body operator generators [6] . Separably measuring the 2-body generators of T yields 
the same syndrome information as the generators of S and although such measurements 
perturb the logical gauge qubits, they commute with Xl and Zl and thus preserve the 
codeword. This remarkable fact is the motivation behind the "zig-zag" network in figure 
[2j which should be more robust to waveguide loss along much of the probing network. 

The fidelity measure used in figure [3] was calculated from p t at each time t as 
Tr(Fp 4 ), where F is the "fidelity operator" constructed by taking the initial state p , 
completely depolarizing all four "gauge" and four relay qubits, and multiplying by 2 4+4 . 
Thus, the expectation of F compares the state at time t to that at 0, tracing over all 
but the logical qubit in the "no-error" syndrome subspace. 
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